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ABSTRACT 

We discuss the implications of rapid (few- minute) variability in the TeV flux of blazars, 
which has been observed recently with the HESS and MAGIC telescopes. The vari- 
ability timescales seen in PKS 2155-304 and Mrk 501 are much shorter than inferred 
light-crossing times at the black hole horizon, suggesting that the variability involves 
enhanced emission in a small region within an outflowing jet. The enhancement could 
be triggered by dissipation in part of the black hole's magnetosphere at the base of 
the outflow, or else by instabilities in the jet itself. By considering the energetics of 
the observed flares, along with the requirement that TeV photons escape without pro- 
ducing pairs, we deduce that the bulk Lorentz factors in the jets must be > 50. The 
distance of the emission region from the central black hole is less well-constrained. We 
discuss possible consequences for multi- wavelength observations. 

Key words: accretion, accretion discs — galaxies: active — BL Lacertae objects: 
individual: PKS 2155-304, Mrk 501 — galaxies: jets — gamma-rays:observations 



1 INTRODUCTION 

The discovery of variable GeV emission from blazars, by 
Compton Gamma-Ray Observatory's EGRET instrument 
(Hartman et al 1992, 2001), opened up a new way to study 
the speed, composition, and energetics of relativistic jets. 
The multi-day variability timescales measured by EGRET 
JL" ' fit comfortably within the prevailing paradigm that vari- 
. £h ' ablity would be imprinted on the scale of the central black 
. hole's horizon, and therefore t var ~ r g /c = GM/c 3 = 1.4mg 
J_i hr, where r g is the gravitational radius and mg = M/10 9 Mq 
is the black hole mass in fiducial units. Disturbances created 
near the black hole could travel outward with a high Lorentz 
factor F (a combination of bulk and pattern speed), before 
radiating energy at a distance > T 2 r g . According to this 
picture, gamma-rays would be produced at ~ (10 2 — 10 4 )r 9 , 
i.e., in a region approaching the scales where most of the 
radio emission is produced. The required Lorentz factors, 
r < 10, are also consistent with the values inferred from 
radio measurements of superluminal motion and brightness 
temperatures. Moreover, placing the GeV emission region 
this far from the black hole greatly alleviates the problem 
of how the photons escape without producing pairs on the 
soft photon background. 

This view has now been challenged by the results of TeV 
observations, which indicate strong variability in at least two 
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blazars (PKS 2155-304: Aharonian et al. 2007; and Mrk 501: 
Albert et al. 2007a), on timescales as short as a few minutes. 
Given the inferred black hole masses of ~ 1O 9 M , these 
timescales are one to two orders of magnitude shorter than 
the shortest timescales expected. Although ultrarelativistic 
motion toward the observer can preserve a short variability 
timescale even when the emission region is far from the black 
hole, it cannot shorten the variability timescale imprinted by 
a source that is stationary in the observer's frame, without 
implausible fine-tuning. Therefore, the TeV results indicate 
that the observed variability is imprinted either by a small 
fraction of the black hole's horizon, or by small-scale fluctu- 
ations intrinsic to the jet itself. 



In this paper we adopt the view that, irrespective of how 
the variability is triggered, it is the jet itself that is producing 
the TeV flares, and study the implications of short-timescale 
variability. In section 2 we present basic scaling relations 
that govern the size and energetics of the flaring regions, 
and we place these constraints in the context of possible 
emission mechanisms in section 3. We show that the large 
apparent luminosity of the flares, and their short timescales, 
constrain the energy content of the emitting regions. These 
constraints, and the requirements that the gamma-rays es- 
cape (section 4), indicate that the flaring regions have bulk 
Lorentz factors > 50, and most likely produce TeV gamma- 
rays via Comptonization of external radiation. 
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2 SIZE AND ENERGETICS OF FLARING 
REGIONS 

We assume that each flaring region, which has a size £' in 
the jet comoving frame, is causally connected during the 
flare (of comoving duration At'), implying £' < cAt' . In the 
lab frame the timescale is dilated to TAt' , but from the ob- 
server's point of view the duration of the flare is compressed 
by a factor (1-/3 cos #) _1 , where j3 = v/c is the dimension- 
less speed and 9 is the angle of motion with respect to the 
line of sight. In the limit /3 « 1 and 6 < T -1 , the compres- 
sion factor approaches 2F 2 . There is a strong observational 
selection effect that favors this limit for the fastest and most 
luminous flares, and we will assume that it holds in the fol- 
lowing analysis. Thus, an observed variability timescale t val 
implies 



var 

c r. 



(i) 



In the lab frame, the total energy content of the flar- 
ing region is related to the comoving energy density e' by 
E ~ e'l'^Y. If a fraction / of this energy is radiated dur- 
ing the flare, the observed power is P b s ~ fE/t var and the 
observer deduces an equivalent isotropic luminosity (i.e., lu- 
minosity inferred from the observed flux without accounting 
for intrinsic anisotropy or beaming effects) of 



4P obs r 2 < 4/e't 2 



c 3 r 6 . 



(2) 



The factor of F 2 in the first relation comes from the fact 
that the power is beamed in the direction of the observer. 
Measured values of Li so can be quite large; for example, 
Li so > 10 46 erg s _1 for the bright TeV flares observed in 
PKS 2155-304 (Aharonian et al. 2007). We will therefore 
normalize Li so to 10 46 erg s _1 . However, the observed TeV 
spectra are quite steep (typical photon spectral indices ~ 
2.5 — 3.5), while EGRET spectra suggest that blazar spectra 
often peak in the 10-100 GeV band. Therefore, if the GeV 
emission flares as rapidly as the TeV emission, Li so could 
well be larger than our fiducial value. 

The shortest variability timescales that have been mea- 
sured to date in PKS 2155-304 (Aharonian et al. 2007) and 
Mrk 501 (Albert et al. 2007a) are 3-5 minutes, so we set 
tvar = SOOts s. We then obtain a lower limit to the internal 
(comoving) energy density in the flaring region, 



i n 9r-lr j.-2t->-6 -3 r-lrr' 

e > 10 / L4et 5 F erg cm ~ / U r , 



(3) 



where U' r is the comoving radiation energy density associ- 
ated with the flare. Note that U' r is smaller than the internal 
energy density. 

We can assess whether the inferred internal energy den- 
sity is reasonable by comparing it to the energy density as- 
sociated with the jet flow, 



(4) 



c r 2 r 2 n 



where Lj is the jet power, r is the distance from the black 
hole and f2 is the solid angle subtended by the jet at r. 
Given the large values of V that we will deduce for the flaring 
regions, we anticipate that jet opening angles may have to be 
be much larger than ~ in order to explain the statistics 
of observable sources. We therefore normalize Q to 0.1 sr. 
Setting x = r/r g , we obtain 



T > 1, 



List 



1/4 



1/4 



1/2 



(5) 



where x& denotes the location of the flare. This relation 
places a significant constraint on F only if the rapid flares 
are produced at radii > 10 3 — 10 4 r s . 

If the flaring regions are indeed moving outward with 
the bulk Lorentz factor of the flow, then 



r 3 + 

£fl > , T = 6 x 10 m Q 



GM 



^r 2 



(6) 



For the values of V estimated below (F > 50), this would 
imply xg > 100. However, it is also possible that the flaring 
regions are patterns that are fixed relative to the lab frame 
(e.g., associated with some external disturbance such as the 
funnel of the accretion flow), in which case the emitting gas 
could be located closer to the black hole. If the flares are 
moving outward, we can combine equations © and ([SJl to 
place a lower limit on the jet power, which is independent 
of the flare timescale: 



Lj > 1.4L iso 



(7) 



3 RADIATION MECHANISMS 

The double humped SEDs of blazars are generally attributed 
to the superposition of a synchrotron spectrum — peak- 
ing in the IR-optical or UV-X-ray for "low-peaked blazars" 
(LBLs) and "high-peaked blazars" (HBLs), respectively — 
and an inverse Compton spectrum (peaking at gamma-ray 
energies) produced by the same population of electrons. 
Both PKS 2155-304 and Mrk 501 are HBLs, as indeed are 
all but one of the blazars detected at TeV energies, to date. 

The seed photons for Comptonization are provided 
primarily by either the synchrotron photons themselves 
(the synchrotron self-Compton [SSC] mechanism: Maraschi, 
Ghisellini & Celotti 1992; Bloom & Marscher 1996) or a ra- 
diation field impinging on the jet from outside (the External 
Radiation Compton [ERC] mechanism: Begelman & Sikora 
1987; Melia & Konigl 1989; Dermer & Schlickeiser 1994; 
Sikora, Begelman & Rees 1994). Since the flare radiation 
density in the comoving frame, U' r (eq. [3]), includes both 
the synchrotron and Compton components, the dominance 
of the Compton (gamma-ray) hump strongly favors the ERC 
mechanism over the SSC mechanism. On the other hand, 
one cannot make such a strong statement if the synchrotron 
peak dominates; in this case, either mechanism is viable. The 
SEDs of the two highly variable TeV sources do not exhibit 
a trend: in Mrk 501 the synchrotron peak appears to dom- 
inate (Albert et al. 2007a), while the Compton component 
dominates the SED of PKS 2155-304 in data presented by 
Foschini et al. (2007; but not in the [non-simultaneous] data 
quoted by Ghisellini et al. 1998). 

The fact that the spectra are quite flat (or even in- 
verted) longward of each peak, and quite steep shortward, 
suggests that most of the energy in the accelerated electrons 
is contained in particles with high random Lorentz factors, 
radiating close to the peak. This view contrasts with ear- 
lier assertions that the peak is associated with cooling of 
the electrons (Sikora et al. 1994; Ghisellini et al. 1998), and 
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could mean that the acceleration mechanism is "particle- 
starved", in the sense that the Poynting flux exceeds the 
kinetic energy flux in the flaring region (Sikora et al. 2005) . 
If particle acceleration is efficient in these flares (i.e., if / is 
not too small and most electrons are accelerated), then we 
may assume that the intrinsic energy density is primarily 
magnetic, e ~ B' 2 /8tt, with 



B' > 2 x 10 6 / -1/2 ^e 2 ^ 1 r _3 G. 



(8) 



According to this emission model, the intrinsic synchrotron 
emissivity of PKS 2155-304 peaks at u Byn ~ W 16 v 16 /r Hz, 
which requires the random Lorentz factors of electrons con- 
tributing to the peak to satisfy 



. n r ft 1/2 rl/4 r -1/4,1/2-p 

7 pcak < 250^ / L 46 t 5 ' T. 



(9) 



If the intrinsic gamma-ray spectrum, peaking at ~ I0 24 /T 
Hz, is produced by Comptonization of the synchrotron spec- 
trum, then 7 P eak ~ 10 4 and F > 50. In this picture, scat- 
tering in the Klein-Nishina regime could contribute to the 
steepness of the TeV spectrum. 

Although large values of V may be needed to produce 
rapidly fluctuating gamma-rays, they can also inhibit effi- 
cient synchrotron cooling of the flare plasma. We have im- 
plicitly assumed that all of the dissipated energy is radiated 
away during the flare. This implies that the cooling timescale 
in the comoving frame is shorter than rt var — if this were 
not the case then our energetic requirements would have to 
increase. If cooling is dominated by synchrotron losses, then 
a sufficient condition to ensure efficient cooling is 



7 > W-*fL^t 5 T 5 



(10) 



To ensure that all electrons with 7 > 7 pca k are able to cool, 
we require 



r < 40^ 1 1 6 /8 r 3/16 ,L^ 1 %- 1/8 - 



(11) 



Thus, an SSC model for rapid flares from PKS 2155-304 can 
barely satisfy the condition for efficient cooling above 7 pea k, 
given our deduction that F > 50; this constraint will become 
tighter once we consider the opacity due to pair production. 

Constraints on radiative efficiency are relaxed consider- 
ably if the gamma rays are produced by the ERC mecha- 
nism. In order for ERC to dominate over SSC, the ambient 
radiation energy density must exceed the synchrotron energy 
density as measured in the comoving frame. The external ra- 
diation density in the lab frame must then satisfy 



[7ext > I0 9 L46^ 2 r~Vg cm" 3 , 



(12) 



if the intensity at xn is approximately isotropic, corre- 
sponding to a luminosity L cxt ~ 2 x I0 35 (r/50) _8 m|a;| 
erg s _1 . Henceforth we normalize P to 50 because both 
energetic and transparency constraints will demand such 
values. If the external radiation illuminates the flaring re- 
gion from behind, subtending a small solid angle f2 ex t (with 
r -2 <C f2 e xt/27T <C 1) in the lab frame, then the required 
energy density and luminosity are increased by a factor 
(fi e xt/27r) -2 . For radiation emitted at r cxt <ra — from the 
central region of an accretion disk, for example — this fac- 
tor is ~ (r cx t/r cxt ) 4 , much steeper than the expected varia- 
tion of emissivity with radiation in most accretion scenarios. 
Therefore, it seems most likely that the external radiation 
responsible for Comptonization is emitted at radii compara- 



ble to Xfi (Dermer, Schlickeiser & Mastichiadis 1992; Dermer 
& Schlickeiser 1993). 

If the (~ isotropic) external radiation field peaks at a 
frequency v cyi t, then the condition for producing the gamma- 
ray peak at v-, ~ 1 24 Hz is 



^oxt 



p2 2 

^ 7pe a k" 



Combining this condition with eq. (]9}, we obtain 



i^ext > 2 x fO 12 ^ 1 / 



-1/2 j 1/2,. 



Hz. 



(13) 



(14) 



Thus, an ERC model for rapidly varying TeV flares would 
require an external radiation source in the submillimeter 
band, if the external radiation is diffuse, and a factor ~ 
(fl ext /2"7r) -1 higher frequency if it comes from behind. 

In contrast to the SSC mechanism, the cooling of rela- 
tivistic electrons in the ERC model becomes more efficient 
with increasing T, with the least efficient cooling occurring 
when the two mechanisms are comparable. For electrons 
near 7 poa k the ratio of cooling time to variability time is sen- 
sitive to r, varying oc T 4 when SSC dominates and oc T~ 4 
when external radiation controls energy loss. 



4 ESCAPE OF RADIATION 

The most stringent conditions on our flare model are set by 
the requirement that the gamma-rays escape without being 
absorbed in 7 — 7 pair production. In blazar models there 
are two possible targets for pair production: the synchrotron 
radiation intrinsic to the jet and the external ambient radi- 
ation. In the case of interactions with jet radiation, optimal 
conditions for escape occur if a TeV photon encounters radi- 
ation only within the flaring region in which it was produced. 
If the TeV photon has to pass through other radiating re- 
gions, this will decrease its escape probability; but given the 
amplitude of observed flares it is plausible that the escape 
constraint is set locallyQ In considering pair production on 
ambient radiation, on the other hand, one must integrate 
over a path length of order r. 

To estimate the pair production opacity internal to the 
flare region, we use the estimate of radiation density U' r from 
eq. ([3} • We assume that this energy density is dominated by 
synchrotron radiation peaking at a frequency 10 16 ^i6/F Hz, 
which therefore has a number density 



fpeakn'„ (fp ea k) ~ 6 x I0 19 i/ 16 1 L46i5 2 r 5 ph s 1 . 



(15) 



The cross section for pair production peaks at <r 77 ~ or/5 
close to threshold, where <jt is the Thomson cross sec- 
tion, and declines at higher energies. Blazar spectra tend 
to be quite flat longward of the synchrotron peak, with 
an energy spectral index a ~ 0.3 — 0.5 (where the flux is 
F v oc v~ a ). Shortward of the peak they decline somewhat 
more steeply than v~ x . Under these circumstances, the most 
probable pair-producing reactions are those close to thresh- 
old. The likely targets for photons with energy I TeV/T 



1 Note, however, that in the Sikora et al. (1994) ERC model for 
3C 279, the pair production constraint due to internal synchrotron 
radiation is much less severe than that due to external diffuse 
radiation. 
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(i.e., the photons with observed energy ~ 1 TeV) have fre- 
quencies ^target ~ 4 x 10 13 F Hz. We therefore need to correct 
fpoakrC^peak) by a factor (f P cak/ftargct) a to account for the 
ratio of target photons to photons near the spectral peak. 
Anticipating that the minimum P for PKS 2155-304 will be 
large enough that target > ^ pC ak, we adopt a = 1. The cor- 
rection factor is then ~ lGOi^ieP -2 . Multiplying by a path 
length £' = ct V!iI F and the threshold cross section, we obtain 
the optical depth to pair production, 

r 77 ~ 2 x 10 10 L 46 ^ 1 r~ 6 . (16) 
The condition t 77 < 1 then implies 

r > 5QL%% 1/6 (17) 

(Celotti et al. 1998). Thus, the escape of TeV photons from 
the site of a rapid flare requires F > 50, provided that the 
observed synchrotron emission comes from the same site. 
This result is independent of the distance of the flaring re- 
gion from the black hole. 

In the presence of diffuse external radiation, the thresh- 
old target frequency for pair production by a 1 TeV photon 
is ~ 6 x 10 13 Hz, or about an order of magnitude higher than 
ZAsxt given by eq. (|14[) , if F ~ 50. Using eq. (|12[1 to estimate 
the external radiation density, and conservatively assuming 
a flat spectral index (a — 0), we estimate a target pho- 
ton number density v syA n v 0>e*t) ~3x io 8 L 46 ^ 2 (r/50)- 8 
ph cm -3 . To obtain the pair production optical depth we 
multiply by the threshold cross section and a path length 
~ r — 1.5 x 10 14 mgx: 

-3 -2 ( r \ ~ 8 
r 77 , tlt ~6x 10 L 4 et 5 I —J mgxa. (18) 

This relation implies that the rapid TeV variability could 
also be produced by Comptonization of diffuse external ra- 
diation, but only out to a radius of several hundred r g for 
r ~ 50. Higher values of F would allow the flare to occur at 
larger radii. 

If the external radiation came from behind the jet, this 
constraint would not be much changed. In a typical pair- 
producing collision, the angle between a TeV photon and 
the soft target would be small, ~ (f2 ox t/27r) 1 ^ 2 , but because 
the incident soft photons must be more energetic by a fac- 
tor (Sl e xt/27r) _1 , the pair production threshold is still about 
an order of magnitude above i^oxt- The ambient radiation 
density required by the ERC model is larger by a factor 
(fl ext /27r) _2 , but this increase is compensated for by the de- 
crease in target photon number and collision rate per target 
photon, each of which scales as fl<, x t/27r. Thus, the optical 
depth for pair production is not much changed, assuming 
that the ERC mechanism produces the TeV flares. If the 
SSC mechanism dominates, so that the external radiation 
supply is weaker, then the pair production constraint would 
be correspondingly relaxed. 

These results are consistent with the ERC model of 
Sikora et al. (1994), which focused on the low-peaked blazar 
3C 279. In that case, the synchrotron peak was at lower 
frequencies, vie ~ 0.1, while the variability timescale was 
taken to be 1 day (ts ~ 300) with F ~ 5. These numbers 
give an external radiation peak frequency in the ultravio- 
let, z^ext ~ 10 16 Hz, and an optical depth for pair produc- 
tion t 77i ex t ~ 7 x 10 _4 m,92:. Thus, flares located at about 
r ~ 10 18 cm are marginally optically thin to pair production. 



5 DISCUSSION 

We have analyzed the requirements for producing rapid 
(few- minute), high amplitude TeV variability in relativistic 
blazar jets. Our model is based on the two standard mod- 
els for gamma-ray emission from blazars, in which the two 
spectral "humps" correspond to synchrotron radiation and 
inverse Compton scattering, respectively. As in the stan- 
dard models, we show that the flaring gamma-rays could 
be produced either by Comptonizing the synchrotron pho- 
tons (SSC model) or by Comptonizing a diffuse background 
source of submillimeter radiation (ERC model). Since the 
flaring regions must be quite small to satisfy causality con- 
straints, they must have large bulk Lorentz factors, F > 50, 
in order for the TeV radiation to avoid pair production 
against the synchrotron photons. VLBI measurements of 
superluminal motion in PKS 2155-304 (Piner & Edwards 
2004) and Mrk 501 (Giroletti et al. 2004) suggest lower val- 
ues of P in these objects, but we stress that the regions 
producing radio emission and gamma-ray flares may have 
very different properties. 

Although both SSC and ERC mechanisms appear capa- 
ble of explaining the rapid TeV variability, ERC seems more 
likely to dominate. Whereas SSC models are tightly con- 
strained (or possibly excluded) between upper limits on F 
imposed by radiative efficiency and lower limits imposed by 
transparency and Comptonization constraints, ERC models 
merely have to satisfy a lower limit on F (since the radiative 
efficiency of ERC increases with F). The radiation energy 
densities required by ERC are modest, and would be hard 
to avoid, even in the presence of a radiatively-inefficient ac- 
cretion flow of the sort likely to be found in BL Lac objects. 
Moreover, the dominance of the Compton peak over the syn- 
chrotron peak, which characterizes the SED of at least one 
of the highly variable TeV sources (PKS 2155-304: Foschini 
et al. 2007), is strong, direct evidence for the ERC process. 

ERC flares could be produced close to the 7 — 7 photo- 
sphere, along the lines suggested by Blandford & Levinson 
(1995). Setting r 77 , cxt ~ 1 in eq. (|18jl and using eq. (|12[1 
to estimate the external radiation density, we obtain an ex- 
ternal luminosity of L cxt ~ 3 x 10 40 LJ 6 2 t 4 (r/50) 8 erg s" 1 . 
This does not represent the entire unbeamed luminosity of 
the blazar, but merely the portion produced at radii > Xfi. 
If the flare occurs at radii ~ 10 2 — 10 3 , as seems likely, 
then such low luminosities could be fully compatible with 
standard accretion models. In particular, the external sub- 
mm radiation could be nonthermal emission produced by 
a radiatively inefficient accretion flow, or thermal emission 
produced by the outer, cool parts of an accretion disc. Al- 
ternatively, it could be produced by relativistic electrons in 
a shear layer surrounding the jet. The effective external lu- 
minosity would vary with radius, depending on the geom- 
etry of the source and the run of F(r). If these behaviors 
were known, one could calculate a relationship between the 
isotropic luminosity of a flare and its duration. The mea- 
sured PDS oc freqr 2 of PKS 2155-304 (Aharonian et al. 
2007) implies Li so oc i var . 

The intrinsic steepness of the TeV spectrum in the 
rapidly varying sources is unlikely to arise from cooling, 
since the spectrum appears to be considerably steeper than 
the hard tail of the synchrotron hump. Klein-Nishina effects 
could contribute to the steepness, as could "self-absorption" 
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due to pair production. The latter effect can occur when 
the TeV photons pair-produce against synchrotron photons 
shortward of the peak. Ignoring Klein-Nishina effects, the 
emissivity for TeV gamma-rays has the same slope as the 
high-energy synchrotron photons (say, j u oc v~ a , where 
q ~ f). But since the absorption coefficient scales with 
TeV frequency as a v oc i/ a , higher energy TeV photons find 
more targets for pair production. Therefore, if the source 
is self-absorbed, the TeV energy spectral index is 2a. This 
would tend to decrease the TeV variability relative to the 
synchrotron photons, if SSC dominates, and might cancel 
the nonlinearity expected from the SSC mechanism. 

Our analysis implies some significant differences be- 
tween LBLs and HBLs that may explain why only the latter 
have exhibited rapid variability at TeV energies (and, with 
the exception of a detection of weak emission from BL Lac 
[Albert et al. 2007b], why only HBLs have been detected in 
the TeV band). Since LBLs have vi6 ~ 10 -4 — 10~ 2 , equa- 
tion ([9| implies that 7 pea k would have to be lower by a factor 
0.01 — 0.1 if these objects were to produce very rapid flares 
via the dissipation of magnetic energy. In this case, gamma- 
rays could not be produced by SSC; ERC would have to 
dominate, with the characteristic frequency of the external 
radiation f ex t in the UV. Such an external radiation source is 
probably not found in BL Lac objects, but might be present 
in optically violently variable (OVV) quasars. In the quasar 
case, however, both the radial scale of the external radiation 
and its luminosity are expected to be quite large, implying 
that the pair photosphere will be located at > 10 4 r 9 . At 
these radii the jet is likely to be so expanded that the energy 
density will be insufficient to create an intense, rapid flare 
(cf. eq. [5]); moreover, the conversion of Poynting flux to ki- 
netic energy (or its loss to radiation) may have largely taken 
place by this radius, further hampering the production of 
flares. Therefore, we might not expect such a system to pro- 
duce very rapid flares with high luminosities. LBL BL Lac 
objects may have trouble producing rapid flares even at GeV 
energies, since they lack a strong source of UV radiation; to 
produce such flares in these objects would require the slow 
dissipation of relatively weak magnetic fields, according to 
equations (|5J and ((9|; thus, we would expect t 5 ^> 1. These 
arguments also suggest a reason for the relative weakness 
of TeV emission, compared to GeV emission, in those LBLs 
that have been detected in gamma-rays. In LBLs, the exter- 
nal Comptonization required to produce TeV photons prob- 
ably extends further into the Klein-Nishina regime, since 
the photon energy in the electron rest frame is roughly the 
geometric mean between the seed photon energy and the 
TeV final energy. This will compound the spectral steepen- 
ing due to effects like self-absorption due to pair production, 
creating a much weaker TeV signal. 

Our model does not make strong predictions about the 
distance from the black hole at which the flares are pro- 
duced. The main constraint is that the flow must have 
accelerated to V > 50 before reaching the flare zone. It 
is well known that relativistic jets, subject to magnetohy- 
drodynamical or fluid forms of pressure, tend to accelerate 
rather slowly. In a ballistic, ultrarelativistic flow F increases 
roughly linearly with radius, suggesting that the jet should 
reach at least 100r 9 before producing the observed flares. If 
the jet is collimated, e.g., by magnetic stresses, the acceler- 
ation could be slower still. However, if the jet is extremely 



"particle-starved" — as we suggested based on evidence that 
the mean random Lorentz factor of electrons in the flare 
zone could be as large as ~ 10 4 — then it may undergo a 
very rapid episode of acceleration close to the black hole. 
By analogy with pulsar winds, magnetocentrifugal stresses 
can accelerate the wind to Lorentz factors ~ cr 1 ^ 3 close to 
the light cylinder, where a is the ratio of energy density to 
rest mass density and the light cylinder would generally be 
expected to lie at a few r g (Michel 1969; Begelman & Li 
1994). If this ratio could be as large as 10 for blazar jets, 
then these jets would probably accelerate to V ~ 100 al- 
most immediately. The remaining energy density would be 
just enough for the electrons to be accelerated to an average 
Lorentz factor of ~ 10 4 . 

The detection of large-amplitude variability on 
timescales shorter than r g /c is important and surprising. We 
can already infer higher Lorentz factors than have generally 
been contemplated for blazars, and that the phenomenon 
is triggered by processes that involve extreme relativistic 
plasmas, perhaps in a black hole magnetosphere. These pos- 
sibilities lend added motivation to future observations - 
and especially to simultaneous multiband observations that 
could discriminate among the options. 
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